Lattice symmetries and regular states in classical frustrated antiferromagnets 
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We consider some classical and frustrated lattice spin models with global 0(3) spin symmetry. 
There is no general analytical method to find a ground-state if the spin dependence of the Hamilto- 
nian is more than quadratic (i.e. beyond the Heisenberg model) and/or if the lattice has more than 
one site per unit cell. To deal with these situations, we introduce a family of variational spin config- 
urations, dubbed "regular states", which respect all the lattice symmetries modulo global 0(3) spin 
transformations (rotations and/or spin flips). The construction of these states is explicited through 
a group theoretical approach, and all the regular states on the square, triangular, honeycomb and 
kagome lattices are listed. Their equal time structure factors and powder-averages are shown for 
comparison with experiments. All the well known Neel states with 2 or 3 sublattices appear amongst 
regular states on various lattices, but the regular states also encompass exotic non-planar states with 
cubic, tetrahedral or cuboctahedral geometry of the T = order parameter. Whatever the details 
of the Hamiltonian (with the same symmetry group), a large fraction of these regular states are 
energetically stationary with respect to small deviations of the spins. In fact these regular states 
appear as exact ground-states in a very large range of parameter space of the simplest models that 
we have been looking at. As examples, we display the variational phase diagrams of the J1-J2-J3 
Heisenberg model on all the previous lattices as well as that of the J1-J2-K ring-exchange model on 
square and triangular lattices. 

PACS numbers: 75.10.Hk,75.40.Cx 



I. INTRODUCTION 

Finding the ground-state (GS) of an antiferromagnetic 
quantum spin model is a notoriously difficult problem. 
Moreover, even classical spin models at zero tempera- 
ture can be non-trivial to solve, unless one carries some 
extensive numerical investigation. In particular there is 
no general method to determine the lowest energy config- 
urations for a simple Heisenberg 0(3) model of the type 

£;-^J(|x,-x,|)S,.S, (1) 

hi 

if the lattice sites {xi} do not form a Bravais lattice. It is 
only if there is a single site per unit cell (Bravais lattice) 
that one can easily construct some GS (see Sec. VI B). 

Another situation where the classical energy minimiza- 
tion is not simple is that of multiple-spin interactions, 
where the energy is not quadratic in the spin compo- 
nents. Finding the GS in presence of interactions of the 
type (Si • Sj)(Sfe • S;) can be difficult and, in general, 
has to be done numerically even on simple lattices with a 
single site per unit cell. Such terms arise in the classical 
limit of ring-exchange interactions. For instance, the - 
apparently simple ~ classical model with Heisenberg in- 
teractions competing with four-spin ring-exchange on the 
triangular lattice is not completely solved." 

In this study, we introduce and construct a family of 
spin configurations, dubbed "regular states". These con- 
figurations are those which respect all the symmetries of 
a given lattice modulo global spin transformations (rota- 
tions and/or spin flips). This property is obeyed by most 



usual Neel states. For instance, a two- (resp. three-) sub- 
lattice Neel state on the square (resp. triangular) lattice 
respects the lattice symmetries provided each symme- 
try operation is "compensated" by the appropriate global 
spin rotation of angle or tt (resp. 0, ±27r/3). 

By definition, the set of regular states only depends 
on the symmetries of the model - the lattice symmetries 
and the spin symmetries - and therefore does not depend 
on the strength of the different interactions (J(|x|) in 
the example of Eq. (1)). These states comprise the well- 
known structures, like the two and three sublattice Neel 
states mentioned above, but also some new states, like 
non-planar structures on the kagome lattice that will be 
discussed in Sec. IV A. 

The reason why these states are interesting for the 
study of frustrated antiferromagnets is that they are good 
"variational candidates" to be the ground-state (GS) of 
many specific models. In fact, rather surprisingly, we 
found that these states (together with spiral states) ex- 
haust all the GS in a large range of parameters of the frus- 
trated spin models we have investigated. For instance, in 
the case of an Heisenberg model on the kagome lattice 
(studied in Sec. VIC) with competing interactions be- 
tween first, second and third neighbors, some non-planar 
spin structures (based on cuboctahedron) turn out to be 
stable phases. In other words, the set of regular states 
and spiral states form a good starting point to determine 
the phase diagram of a classical 0(3) model, without 
having to resort to lengthy numerical minimizations. 
In several cases, we even observed that one of the reg- 
ular states reaches an exact energy lower bound, there- 
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fore proving that it is one (maybe not unique) GS of the 
model. 

These states may also be used when analyzing exper- 
imental data on magnetic compounds where the lattice 
structure is known, but where the values (and range) 
of the magnetic interactions are not known. In such a 
case, the (equal time) magnetic correlations - measured 
by neutron scattering - can directly be compared to those 
of the regular states. If these correlations match those of 
one regular state, this may be used, in turn, to get some 
information about the couplings. With this application 
in mind, we provide the magnetic structure factors of all 
the regular states we construct and powder-averages of 
some of them (see App. A). 

The organization of the paper is as follows: in Sec. II 
we present the definition of a regular state, a state that 
weakly breaks the lattice symmetry and all the nota- 
tions needed for the group theoretical approach. In 
Sec. Ill A, we explain the algebraic structure of the group 
of joint space- and spin-transformations that leave a reg- 
ular spin configuration invariant (Algebraic symmetry 
group) and then explain how to construct regular states 
(Sec. Ill B). This approach is algebraically very similar to 
Wen's construction of symmetric spin liquids, but there 
are also strong differences in the invariance requirements 
(see App. B): whereas the symmetric spin liquids do not 
break lattice symmetries (they are "liquids") , our regular 
states indeed break lattice symmetries but in a "weak" 
way. These sub-sections are self-contained, but can be 
skipped by readers interested essentially in the results. 
In Sec. IIIC, we give an example of such a contruction 
on the triangular lattice and list the regular states on this 
lattice. In sections Sec. IV A and IV B we list the regu- 
lar states on the kagome and honeycomb lattices (which 
have the same algebraic symmetry group as the triangu- 
lar lattice), and with a minimum of algebra we present 
the regular states on the square lattice (Sec. IV C). We 
then show that spiral states can be seen in this picture 
as regular states with a lattice symmetry group reduced 
to the translation group (Sec. IV D). In Sec. V we discuss 
geometrical properties of regular states and the relation- 
ship between regular states and representations of the 
lattice symmetry group. This section can be skipped by 
readers more interested in physics than in geometry. In 
Sec. VI we study the energetics of these regular states 
and therefore their interest for the variational descrip- 
tion of the T = phase diagrams of frustrated spin mod- 
els. We first show in Sec. VIA that all regular states 
which do not belong to a continuous family are energet- 
ically stationary with respects to small spin deviations 
and thus good GS candidates for a large family of Hamil- 
tonians. After having given a lower bound on the energy 
of Heisenberg models (Sec. VI B), we then show that over 
a large range of coupling constants the regular states are 
indeed exact GSs of the J1-J2-J3 model on the honey- 
comb and kagome lattices (Sec. VIC). We then display 
in Sec. VI D a variational phase diagram of the J1-J2-K 
model on square and triangular lattices. In Sec. VIE we 



discuss finite temperature phase transitions: the non pla- 
nar states are chiral and should give rise to a T ^ phase 
transtion. Sec. VII is our conclusion. Powder-averages of 
the structure factors of the regular states on triangular 
and kagome lattices are displayed in App. A. Analogies 
and differences between the present analysis and Wen's 
analysis of quantum spin models are explained in App. B. 

II. NOTATIONS AND DEFINITIONS 

We will mostly concentrate on Heisenberg-like models 
where on each lattice site i, the spin Si is a three compo- 
nent unit vector. But the concept of regular state can be 
easily extended to the general situation where Si belongs 
to an other manifold A (as for example for nematic or 
quadrupolar order parameters). 

We note by Ss the group of the "global spin symme- 
tries" of the Hamiltonian. In the general framework, an 
element of Ss is a mapping of A onto itself which does 
not change the energy of the spin configurations. For an 
Heisenberg model without applied magnetic field, Ss is 
simply (isomorphic to) the orthogonal group 0(3). In a 
similar way, we note by Sl the lattice symmetry group of 
the Hamiltonian. An element of Sl acts on spin config- 
urations by mapping the lattice L onto itself and is the 
identity in the spin space A. 

In this paper, we will restrict ourselves to the (rather 
common) situation where the full symmetry group Sh of 
the model is the direct product Ss Sl- 

Let be the set of all the applications from the lattice 
symmetry group Sl to the spin symmetry group Ss- An 
element G of Q associates a spin symmetry Gx to each 
lattice symmetry X: 

G : Sl ^ Ss /r,^ 
X^Gx ^ ' 

We now concentrate on a fixed spin configuration c. We 
note He its stabilizer, that is the subgroup of Sh which 
elements do not modify c. Its spin symmetry group is 
the group of unbroken spin symmetries: = Ss r\ He- 

Definitions: 

• A mapping G £ Q is said to he compatible with 
a spin configuration c if the composition of an ele- 
ment of Sl with its image by G leaves c unchanged: 

WXeSl, GxXeH,. (3) 

• A configuration c is said to be regular if any lat- 
tice symmetry X € Sl can be "compensated" by an 
appropriate spin symmetry Gx G Ss, which means 
GxX\c) ~ |c) (that is GxX E He). In other words, 
c is regular if there exists a mapping G € G such 
that G and c are compatible. 

In a regular state, the observables which are invariant 
under Ss are therefore invariant under all lattice sym- 
metries. These definitions are summarized in Fig. 1. 



FIG. 1: (Color online) A lattice symmetry X £ Sl acts on a 
spin configuration c to give a new configuration c' = Xc. If c 
is regular, there is a spin symmetry Gx G Ss such that one 
gets back the initial state: Gxc' = c. 

The simplest regular states are those which are already 
invariant under lattice symmetries (i.e. Sl C H^), with- 
out the need to perform any spin symmetry. This is the 
case of a ferromagnetic (F) configuration, with all spins 
oriented in the same way. But less trivial possibilities 
exist, as the classical GS of the antiferromagnetic (AF) 
first neighbor Heisenberg interaction on the square lat- 
tice. This GS possesses two sublattices with opposite spin 
orientations. Each lattice symmetry X either conserves 
the spin orientations, or reverses them, so we can choose 
as Gx either the identity or the spin inversion — > — S^. 

If the subgroup = Ss He of unbroken spin sym- 
metries contains more than the identity, there are several 
elements of Q compatible with c. For each X, they are 
as many Gx as elements in . In the previous exam- 
ple of the GS of the AF square lattice, is the set of 
spin transformations that preserve the two opposite spins 
orientations: this group is isomorph to 0{2). Beginning 
with a compatible G, each Gx can be composed with an 
element of to give an other compatible element of Q. 

To summarize, regular states are not restricted to 
states strictly respecting the lattice symmetries, but to 
states that in some way weakly respect them. We will 
now explain how to construct all the regular spin config- 
urations on a given lattice. 



III. CONSTRUCTION OF REGULAR STATES 

To construct the regular states, we proceed in two 
steps. In the first step, we fix a given unbroken spin sym- 
metry group H^, and consider the algebraic constraints 
that the lattice symmetry group Sl imposes on a map- 
ping G € Q, assuming that some (so far unknown) spin 
configuration c is compatible with G. These constraints 
lead to a selection of a subset Q"^ of G, composed of the 
mappings G which are compatible with the lattice sym- 
metries. For an element G of the group 

H° = {GxX, XeSL}x H^, (4) 

is dubbed the algebraic symmetry group associated to G. 

In the second step, one determines the configurations 
(if any) which are compatible with a given algebraic sym- 
metry group. 
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FIG. 2: (Colour online) Generators of the lattice symmetries 
group Sh for the triangular, kagome, honeycomb and square 
lattices. For the first three lattices : the two translations T\ 
and T2 (along the two basis vectors Ti and T2), the reflexion 
a an the rotation R(, of angle 7r/3. For the square lattice, 
generators of 5*^ are Ti, T2-, o an the rotation _R4 of angle 
7r/2. 



A. Algebraic symmetry groups 

We fix the spin symmetry group (to be exhaus- 
tive, we will consecutively consider each possible i?^). 
Let X, Y and Z, three elements of Sl such that XY — 
Z. We will see that this algebraic relation imposes 
some constraints on the mappings G which are com- 
patible with a spin configuration. Indeed, we assume 
that there exists a configuration c compatible with G. 
Then, GzZ and GxXGyY are in He- This implies that 
GxXGyYZ-^Gz^ is also in He- Elements of Sl and Ss 
commute, so we have GxXGyYZ^^G^^ — GxGyG^^ , 
which is a pure spin transformation. We deduce that 

•iX,YeSL, GxGyGl^y^eH^'. (5) 

If the constraint above is not satisfied, G must be ex- 
cluded from the set Q"^ of the algebraically compati- 
ble mappings. contains only elements of Q verifying 
Eq. (5). 

Now we illustrate these general considerations using 
the following example: L is an infinite triangular lattice 
and the spin space A is the two-dimensional sphere S2 
(Heisenberg spins). Sl is generated by two translations 
Ti and T2 along vectors Ti and T2, a reflexion a and a 
rotation Rq of angle 7r/3, described in Fig. 2 and defined 
in the (Ti,T2) basis as: 

Ti : (ri,r2)^(ri + l,r2) (6a) 

T2 : (ri,r2)^(ri,r2 + l) (6b) 

a : (ri,r2) (r2,ri) (6c) 

R6-{ri,r2)'-^{ri~r2,ri). (6d) 

The spin symmetry group Ss is chosen to be 0(3) (as 
for an Heisenberg model). In such a system, the unbro- 
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ken symmetry group is either isomorph to {/}. Z2 
or 0(2), depending on the orientations of the spins (non- 
coplanar, coplanar or colinear respectively). The non- 
planar case, = {/}, is the most interesting case and 
we choose it for this example. The two other cases can 
be treated by reducing A to the circle Si or 5o = {1,-1} 
{XY or Ising spins) and 5*5 to 0{2) or 0(1) in order to 
have = {I}, which considerably simplifies the calcu- 
lations. 

We assume that a mapping G belongs to (alge- 
braically compatible). As — {/}, Eq. (5) allows to 
construct the full mapping G simply from the images of 
the generators of the lattice symmetry group Sl- As 
several combinations of generators can produce the same 
element of Sl, the images by G of the Sl generators 
must satisfy some algebraic relations. These relations 
where needed in a similar algebraic study in Ref. and 
consist in all the relations necessary to put each product 
of generators in the form a" RqT^^ , where s = 0, 1, 
r = 0, 1, 5 and ti, t2 G Z. These relations are: 



TiRqT2 — i?6 

R6T1T2 = T2Re 
Tier = aT2 

RquRq — a. 

From these equations and from Eq. (5) we get: 



(7a) 
(7b) 
(7c) 
(7d) 
(7e) 
(7f) 
(7g) 



are listed below : 

Gti = ~ l,Ga — Sal, G Be ~ SrI, 

Gti — Gt2 — Ii Ga — £a-I, Glie — ^nfizTT, 
Gti = Gt2 ~ I, Ga = £aRzTr, G Rg = £rI, 
Gti = Gt2 — Ii Gc = SaRzTT, Gr^ — SrRztt, 
Gti = Gt2 = I, Ga = SaRziT, G Rg ~ SrRxtt, 
Gti — Gt2 = Ii Ga = SaRzTT, G Rg ~ SrRxO, 

010^ 

RxTT, Gt2 = RyiT, Ga = —Sa 100 



.001, 



Gi 




Gti — Gt2 — Rzi^i Ga ~ £crl, Grq — eRRx.Tr, 
Gti = Gt2 = -Rz4f J Gcr — £a-RziT, Gr^ — SrRxtt- 



(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 
(lOf) 



(lOg) 



(lOh) 
(lOi) 



where x _L z, Sa, £r G {—1, 1} and 9 £ {f 7 if }• Each 
line corresponds to 4 solutions, except Eq. (lOf) with 8 
solutions. We stress that the algebraic symmetry groups 
depend on Ss, He and on the algebraic properties of Sl, 
but not directly on the lattice L. In particular, different 
lattices can have the same algebraic symmetry groups. 
The results Eqs. 10 are exactly the same on a honeycomb 
or a kagome lattice with symmetries of Fig. 2 because the 
algebraic equations Eqs. 7 stay the same. 



Gti Gt2 


— Gt2 Gti 


(8a) 


Gti Grq Gt2 


= GRe 


(8b) 


Grq Gti Gt2 


= Gt2 GRe 


(8c) 


GtiGc 


— GcrGja 


(8d) 


'-'Re 


= / 


(8e) 


Ga^ 


= / 


(8f) 


G ReGaG Rg 


— G(j. 


(8g) 



To solve this system, we first remark from Eq. (8a) and 
(8d) that Gti and Gtj commute and can be obtained 
from each other by a similarity transformation. Thus, 
we are in one of these four cases 



Gti — Gt2 — I, 
9ti = Ot2 = I" and -L , 

Gti = Gt2 7^ I, 
Gt, ^G^j,^ and Ot^ ^ tt. 



(9a) 
(9b) 
(9c) 
(9d) 



There, we have labeled the elements of 0(3) by a rotation 
of axis n and an angle 6 G [0, tt] (times a determinant 
e = ±1, not appearing here). Up to a global similarity 
relation, we obtain 28 solutions of the system of Eqs. (8) 
in the case of Eq. (9a), 4 for Eq. (9b), 8 for Eq. (9c) and 
for Eq. (9d) (it contradicts Eq. (8c)). The 40 solutions 



B. Compatible states 

The second step consists in taking each element of 
and finding all the compatible states. This last step is 
fully lattice dependent. 

To construct a regular state compatible with some 
mapping G G Q^, one first chooses the direction of the 
spin on a site i. Then, by applying all the transforma- 
tions of Sl, we deduce the spin directions on the other 
sites. A constraint appears when two different transfor- 
mations X and Y lead to the same site X{i) — Y{i). The 
image spins have to be the same: Gx{Si) = Gy(Si). It 
can either give a constraint on the direction of Sj , either 
indicate that no G— compatible state exists. 

To find these constraints, we divide the lattice sites in 
orbits under the action of Sl (if all the sites are equiv- 
alent, there is a single orbit). In each orbit, we choose 
a site i. Each non trivial transformation X that does 
not displace i gives a constraint: Gx(Sj) = Si. For each 
G G G^, the associated regular states are obtained by 
choosing a site in each orbit, a spin direction respecting 
the site constraints and then propagating the spin direc- 
tions through the lattice using the symmetries in Sl- 



C. Example of regular state construction: the 
triangular lattice 

Let us apply this method to the example of the triangu- 
lar lattice. There is a single orbit, and only two distinct 
and non trivial transformations leave invariant the site 
of coordinates (0,0) in the (Ti,T2) basis (see Fig. 2): 
a and i?6" for 1 < n < 5, giving the two constraints 
G,(S(0,0)) = Gj^,(S(0,0)) = S(0,0). 

The mapping of Eq. (10a) has compatible states only 
for Eli — £(7 = ^- They are ferromagnetic (F) states, as 
shown in Fig. 3(a). Since Gti_2 = ^ for the Eqs. (10b)- 
(lOf), no new regular states can be compatible with any 
of them. 

The mapping of Eq. (lOg) has compatible states only 
for £r = 1 and = -1. Then S(0,0) = ±(1,1,1)/%/! 
and the state is the tetrahedron state depicted in 
Fig. 3(b), where the spins of four sublattices point to- 
ward the corners of a tetrahedron. The sign of S(0, 0) 
determines the chirality of the configuration. 

The next regular state is the coplanar state of Fig. 3(c), 
which is compatible with Eq. (lOi) for — Sa = ^ and 
S(G,0) = ±(1,0,0). The three sublattices are coplanar 
with relative angles of 120°. This state is not chiral be- 
cause the configurations obtained with the two possible 
S(0,0) are related by a global spin rotation in SO{2>). 

A continuum of umbrella states are compatible with 
Eq. (lOi) with e„ = 1 and eji — —1. They are depicted 
in Fig. 3(d), where the sublattices are the same than for 
the coplanar states but the relative angles between the 
spin orientations are all identical and < 120°. This family 
interpolates between the F and the coplanar states. 

We started by choosing = {/}, but states with 
ijS ^ 22 (for the coplanar state) or 0(2) (for the F 
state) have been obtained anyway. One can check that 
choosing another would not give any new regular 
state. All the regular states are thus those gathered in 
Fig. 3. 

The Bragg peaks of these states are displayed in the 
hexagonal Brillouin zone in the right column of Fig. 3 
and their powder-averaged structure factors in App. A 
together with the formulas for these quantities. 



IV. REGULAR STATES FOR HEISENBERG 
SPINS ON SEVERAL SIMPLE LATTICES 

In the following we enumerate the regular states on 
the kagome and honeycomb lattices, two lattices which 
have a symmetry group Sl isomorphic to that of the 
triangular lattice. To be complete, we also present the 
regular states on the square lattice and discuss the spiral 
states that may be seen as regular states when reduces 
to the translation group. 
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(a) Ferromagnetic (F) state. E = 6Ji + 6J2 -I- 6J3 -I- i2K. 




(b) Tetrahedral state. E = -2Ji - 2 J2 + 6J3 - 3iK/3. 




(c) Coplanar state. E = — 3Ji -I- 6J2 — SJs — SK. 




(d) F umbrella states 



FIG. 3: (Color online) Regular states on the triangular lattice. 
The sublattice arrangements (labelled by colors) and the spin 
directions on each sublattice are displayed in the left and cen- 
ter columns. A spin unit cell is surrounded with green lines. 
The positions and weights of the Bragg peaks in the hexago- 
nal Brillouin zone of the lattice are in the right column. The 
energy per site of each structure is given as a function of the 
parameters of the models described in Sec VI. 



A. Kagome lattice 

The symmetry group Sl of the kagome lattice is iso- 
morphic to that of the triangular lattice, thus the alge- 
braic solutions Eqs. (10) remain valid. Carrying out the 
approach of Sec. IIIB for this new lattice, one obtains 
all the regular states on the kagome lattice. They are 
displayed in Fig. 4 with the positions and weights of the 
Bragg peaks. The equal time structure factor is depicted 
in the Extended Brillouin Zone (EBZ), drawn with thin 
lines in Fig. 4: the kagome lattice has 3 spins per unit 
cell of the underlying triangular lattice and the EBZ has 
a surface four times larger than the BZ of the under- 
lying triangular Bravais lattice, drawn with dark lines. 
Powder-averaged structure factors of the regular states 
are given in App. A. 
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(a) Ferromagnetic (F) state. 
£; = 4Ji+4J2+24 + 4J3. 




(b) q = state. E = -2Ji - 2J2 + 2J^ + 4Ja. 




(d) Octahedral state. E = 2J!^ - 4Ja. 




(e) Cubocl state. E = -2Ji + 2J2 - 2J^. 



One regular state is colinear (fff = 0(2)): 

• the ferromagnetic (F) state of Fig. 4(a). 

Two states with a zero total magnetization are copla- 
nar (fff = Z2): 

• the q = state of Fig. 4(b) has 3 sublattices of 
spins at 120° and a 3 sites unit cell, 

• the \/3 X y/3 state of Fig. 4(c) has 3 sublattices of 
spins at 120° and a 9 sites unit cell. 

Three states with a zero total magnetization com- 
pletely break 0(3) (iJf = {/}): 

• the octahedral state of Fig. 4(d) has 6 sublattices of 
spins oriented toward the corners of an octahedra 
and a 12 sites unit cell, 

• the cubocl state of Fig. 4(c) has 12 sublattices of 
spins oriented toward the corners of a cuboctahe- 
dron and a 12 sites unit cell, 

• the cuboc2 state of Fig. 4(f) has 12 sublattices of 
spins oriented toward the corners of an cuboctahe- 
dron and a 12 sites unit cell. Note that the first 
neighbor spins have relative angles of 60°, in con- 
trast to 120° for the cubocl state. 

Two continua of states with a non-zero total magneti- 
zation completely break 0(3) {H^ = {/}): 

• the q = umbrella states of Fig. 4(g), left, 

• the \/3 X \/3 umbrella states of Fig. 4(g), right. 

These continua interpolate between the ferromagnetic 
state and the coplanar states Fig. 4(b) and Fig. 4(c). 



B. Honeycomb lattice 




< r > 

1/6 1/6 



(f) Cuboc2 state. E = 2Ji- 2J2 - 2^. 





q = (left) and \/3 X \/3 (right) umbrella states 



FIG. 4: (Color online) Regular states on the kagome lattice 
and their equal time structure factors in the EBZ (see text). 
The energies (per site) of these states are given for the J1-J2- 
J3-J3 model described in Sec VI. 



All the regular states on the honeycomb lattice are 
depicted in Fig. 5 and listed below. The EBZ is drawn 
with thin lines (its surface is three times larger than that 
of the BZ). 

Two regular states are colinear (ijf = 0(2)): 

• the ferromagnetic state of Fig. 5(a), 

• the antiferromagnetic state of Fig. 5(b) has 2 sub- 
lattices of spins oriented in opposite directions and 
a 2 sites unit cell. 

Two states with a zero total magnetization completely 
break 0(3) (i/f = {/}): 

• the cubic state of Fig. 5(c) has 8 sublattices of spins 
oriented toward the corners of a cube and a 8 sites 
unit cell, 

• the tetrahedral state of Fig. 5(d) has 4 sublattices of 
spins oriented toward the corners of a tetrahedron 
and a 4 sites unit cell. 
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(a) Ferromagnetic (F) state. E = 3Ji + 6J2 + 3J3. 




(b) Antiferromagnetic (AF) state. E = — 3Ji + 6J2 — 3J3. 




(c) Cubic state. E = Ji - 2J2 - 3J3. 




(d) Tetrahedral state. E = -Ji - 2J2 + 3J3. 




(e) V states 



FIG. 5: (Color online) Regular states on the honeycomb lat- 
tice and their equal time structure factors in the EBZ (see 
text). The energies (per site) are given for a J1-J2-J3 Heisen- 
berg model. 



A continuum of states with a non-zero total magneti- 
zation partially breaks 0(3) {H^ = Z2): 



• the V states of Fig. 5(e), which interpolate between 
the F and AF states. 



C. Square lattice 

The symmetry group Sl of the square lattice is distinct 
from that of the triangular lattice (see Fig. 2) and one 
has to determine its algebraic symmetry groups. They 
are listed below: 

GTi = = Ga = EaRzTTi G = E rRxti, 

Gti = Gt2 = £lRzTrSi,Gcr = Sa-RxiriGji^ — E fiRy,2L , 
Gti = Gt2 — £lRz7rSi,Ga = E a RzirS, , G — ErRztvSr, 
Gti — Gt2 = SlRzTT, Ga = Ea-RzTrS„, G R^ = E rR^tt , 
Gti — Gt2 — SlRzTT, Ga = EaRxn, Gr^ — ErR^ttSr, 
Gti — Gt2 = £lRzTT, Ga ~ Ea-RxTT, G R^ — ErRxtt, 

GTi = Gt2 = £\RzTT^ Gc — E^RxTT, G R^ — £fli?y7r, 

Gti — El^xTT, Gt2 — £lRyin 
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where x, y, z are orthogonal vectors, ei, 62, 63, £1, e^, 
Er — ±1 and Sr, 5^, 5i = or 1. 

Then, the construction of the compatible states leads 
to the regular states depicted in Fig. 6 and listed below. 

Two regular states are colinear (fff = 0(2)): 

• the ferromagnetic state of Fig. 6(a), 

• the (tt, vr) Neel (AF) state of Fig. 6(b) has 2 sub- 
lattices of spins oriented in opposite directions and 
a 2 sites unit cell. 

One state with a zero total magnetization is coplanar 
(Hf =Z2): 

• the orthogonal coplanar state of Fig. 6(c) has 4 sub- 
lattices of spins with angles of 90° and a 4 sites unit 
cell. 

Then we have three continua of states with different spin 
symmetry group : 

• the V states of Fig. 6(d) have a non-zero total mag- 
netization and partially break 0(3) {H^ = Z2). 
They interpolate between the F and the (tt, tt) Neel 
states, 

• the tetrahedral umbrella states of Fig. 6(e) have 
a zero total magnetization and completely break 
0(3) (iJf = {/}). They interpolate between the 
(tt, tt) Neel and the orthogonal coplanar state, 

• the 4-sublattice umbrella states of Fig. 6(f) have a 
non-zero total magnetization and completely break 
0(3) (iJf {/}). They interpolate between the F 
and the orthogonal coplanar state. 
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(a) Ferromagnetic (F) state. 
£ = 4 Ji + 4 J2 + 4J3 + UK. 




(b) (tt, n) Neel (AF) state. 
E = -4Ji + 4J2 + 4 J3 - 2K. 




(e) Tetrahedral umbrella (f) Umbrella states 

states {AF umbrellas) {F umbrellas) 

FIG. 6: (Color online) Regular states on the square lattice 
and their equal time structure factors in the square BZ. The 
energy per site of each structure is given as a function of the 
parameters of the models described in Sec VI. 

D. Regular states with only translations 

When the lattice symmetry group is commutative, the 
construction of regular states is particularly simple. This 



occurs if only translations are considered. In that case, 
one may chose some arbitrary directions for the spins of 
the reference unit cell. Then, one has to chose an 0(3) 
element Gt^ associated to each unit lattice translation 
Ti in direction i (with as many generators as space di- 
mensions). Assuming that — I, and using the fact 
that the translations commute with each other, we find 
that the Gt^ also commute. A first family of solutions 
consists in choosing a set of rotations with the same axis 
n, and unconstrained angles. This gives the conventional 
spiral states. Thanks to the arbitrary choice of the spin 
directions in the reference unit cell, such states are not 
necessarily planar. 

Finally, an other family of solutions can be obtained 
by choosing the Gt^ among the set of 7r-rotations with 
respect to some orthogonal spin directions, therefore in- 
suring the commutativity. 

All these solutions may be generalized by combining 
one or more Gt; with a spin inversion — /. These gener- 
alized spiral states will be noted SS in the following. 

V. GEOMETRICAL REMARKS 

In this section, we discuss some geometrical properties 
of regular states. 

A. Groups and polyhedra 

From a regular state c, one can consider the set 'S C A 
of all the different orientations taken by the spins. We 
assume that c has a finite number of sublattices/spin di- 
rections, so that E is finite. For a three component spin 
system, S is just a set of points on the unit sphere S2, 
as displayed in Figs. 3, 4, 5 and 6. E may be a single 
site, the ends of a segment, the corners of a polygon or a 
polyhedra. 

The four lattices studied here share some special prop- 
erties: all the sites and all first-neighbors bonds are 
equivalent (linked by a 5*^ transformation). Due to this 
equivalence, E also form a segment /polygon/polyhedron 
with equivalent vertices and bonds. A polyhedron with 
this property is said to be quasi-regular. If the elemen- 
tary plaquettes of the lattices are also equivalent (as in 
the triangular, square and hexagonal lattices, but not in 
the kagome lattice where both triangular and hexagonal 
elementary plaquettes are present) and if S is a polyhe- 
dron, its faces should also be equivalent. E must then 
be one of the five regular convex polyhedra (Platonic 
solids):^' tetrahedron, cube, octahedron, dodecahedron 
or icosahedron. ' ' 

We now only consider the case where — {1} (this 
condition can always be verified by reducing A to its el- 
ements invariant by and by consequently modifying 
Ss)- Clearly, the lattice symmetries constraint the possi- 
bilities for the set E, since each lattice symmetry X per- 
mutes the sites in E but leaves it globally unchanged. 



9 



But since the state c is regular, these permutations can 
also be achieved by a spin symmetry in Ss, and the sym- 
metry group S's of E should be viewed as a finite sub- 
group of Ss- 

For 5*5 = 0(3), the classification of these subgroups - 
called point groups - is a classical result in geometry, 
it contains seven groups (related to the three symmetry 
groups of the five regular polyhedra) and seven infinite 
series (conventionally noted C„, C'nv, Cnh, Dn Dnh D„d 
and Sn with n G N. They are related to the cyclic and di- 
hedral groups). Of course, the non planar regular states 
we have discussed so far (Sec. Ill C and IV) fall into this 
classification. For instance, the three- and four- sublat- 
tice umbrella states of Fig. 3(d), 6(e) and 6(f) correspond 
to Csi,, I?2d and (with respectively 6, 8 and 8 ele- 
ments). The cubic, octahedral and cuboctahedronl states 
correspond to the symmetry group of the cube (48 ele- 
ments), and the tetrahedral state corresponds (of course) 
to its own symmetry group. 



B. Regular states and representation of the lattice 
symmetry group 

We again focus on three-component spin systems with 
a spin symmetry group Ss = 0(3). In a regular state c 
each lattice symmetry X can be associated to a matrix 
Gx in 0(3). Now, as in Sec. Ill A, we can compare the 
actions of two lattice symmetries X and Y. CxGyG^y 
belongs to . By an appropriate choice of Gx, it is 
possible to obtain GxGy — Gxy, which implies that G 
is a representation of the lattice symmetry group Sl. Its 
dimension is 1 for a colinear state, 2 for planar states, 
and 3 for the others. Is this representation reducible 
? If yes, it must contain at least one representation of 
dimension 1 (because the maximal dimension considered 
here is 3), thus there exist at least one spin direction 
which is stable under all the spin symmetry operations 
spanned by Gx with X d Sl- Except in the trivial 
colinear case, one can easily check that it is the case 
only for the states belonging to a continuum. For the 
V-states, G is the tensor product of a trivial and a non 
trivial Id representation of Sl (and of any Id third one). 
For the umbrella states, G is the tensor product of a 
trivial Id and a 2d irreducible representation (IR). For 
the tetrahedral state of Fig. 6(c), G is the tensor product 
of a non-trivial Id and a 2d IR representation. For the 
other cases, the associated representation is irreducible. 

There is another context where antiferromagnetic Neel 
states are known to be related to irreducible representa- 
tions. If a quantum antiferromagnet has a GS with long- 
range Neel order, its spectrum displays a special struc- 
ture, called "tower of states". It reflects the fact that 
a symmetry breaking Neel state is a linear combination 
of specific eigenstates with different quantum numbers 
describing the spatial symmetry breaking, and with dif- 
ferent values of the total spin S, describing the SU{2) 
symmetry breaking. If such a quantum system has a 



GS with a regular Neel order, its tower of state should 
have an S' = 1 state with the same quantum numbers 
as those of the irreducible representation X i-> Gx dis- 
cussed above. The reason why this representation shows 
up in the S = I sector of the tower of state is because 
5 = 1 corresponds to the action of the lattice symmetries 
onto a three-dimensional vector, as the classical spin di- 
rections. 



VI. ENERGETICS 

As discussed in the introduction, there is no simple way 
to find the GS of a classical spin model if the lattice is 
not a Bravais lattice, and/or if spin-spin interactions are 
not simply quadratic in the spin components. So far, we 
have discussed regular states from pure symmetry con- 
siderations, but in Sec. VIA we show that, under some 
rather general conditions, a regular state is a stationary 
point for the energy, whatever the Hamiltonian (provided 
it commutes with the lattice symmetries) . 

In addition, we argue that regular states are good can- 
didates to be global energy minima. To justify this, we 
first discuss a rigorous energy lower bound (Sec. VI B) for 
Heisenberg like Hamiltonians and investigate in Sec. VI C 
several Heisenberg models with further neighbor interac- 
tions (Ji, J2, J3, etc.) on non-Bravais lattices such as 
the hexagonal and kagome lattices. In large regions of 
the phase diagrams, one regular state energy reaches the 
lower bound and is one (may be not unique) exact GS. 



A. A condition for a regular state to be 
"stationary" with respect to small spin deviations 

To address the question of energetic stability of regular 
states, we give some conditions under which an infinites- 
imal variation of the spin directions would not change 
the energy (necessary condition to have a GS). To sim- 
plify the notations we consider an Heisenberg model with 
some competing interactions (such as in Eq.( 1)), but the 
arguments easily generalize to multi-spin interactions of 
the form (Si • Sj)(Sfc • S/) . . . (respecting the lattice sym- 
metries) . 

We assume that there is a non trivial lattice symmetry 
X which leaves one site i unchanged: X{i) = i (existence 
of a non-trivial point group). In addition, we assume 
that a spin rotation Rg of axis n and angle 9^0 can 
be associated to X in order to have RgXc — c. These 
conditions insure that the invariant direction of Kg is 
n = ±Si. Excepted states belonging to a continuum, all 
regular states verify these conditions on the lattices we 
have studied . 

With these conditions, the derivatives of the energy 
with respect to the spin directions vanish. The proof is 
as follows. One considers the local field = which 
is experienced by the spin i. is a linear combination 
of the Sj where j runs over the sites which interact with 
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the site i: 

h» = E^'^ E S,' (12) 

where Nd{i) is the set of the neighbors of i at distance 
d on the lattice. Since the configuration c is invariant 
under RgX, one may also compute as 

X reshuffles the neighbors of i (at any fixed distance) 
but since X{i) — i, Nd{i) is globally stable: Nd{i) = 
X{Ndii)). So, from Eq. (13), we have 

h, =i?,(h,). (14) 

We therefore conclude that h,j is colinear with n and thus 
colinear with Si. This shows that the energy derivative 
vanishes for spin variations orthogonal to (longi- 
tudinal spin variations are not allowed as (Sj)^ must be 
kept fixed). 

All regular states studied in the previous examples that 
do not belong to a continuum are thus energetically sta- 
tionnary with respect to small spin deviations. They are 
thus interesting candidates for global energy minima. 

B. Lower bound on the energy of Heisenberg 
models 

The Fourier transform Sq^ of the local spin on a peri- 
odic lattice of N unit cells is defined by 

^1* = T^E^-^""'"- 

where each site is labeled by an index i = 1 ... to (to is 
the number of sites per unit cell), x is the position of 
its unit cell, and q is a wave vector in the first Brillouin 
zone. For an Hamiltonian in the form of Eq. (1), the 
energy can be written as: 

(15) 

V,X 

z J — l,...,m 

= E ^..(q)S-q,-Sq, , (16) 

qeBZ 
— l,...,m 

with 

V 

Since (S,x)^ = 1 for all i and x, Xl^x = E^q S^x = 
mN, we see that a lower bound on the energy (per site) 
is obtained from the lowest eigenvalue of the matrices 
J(q):^ 

^ > min (J™") (18) 
mN ~ {q} ^ ^ 



where J™™ is the lowest eigenvalue of the matrix J(q). 

If the lattice has a single site per unit cell (m = 1) this 
lower bound is reached by a planar spiral of the form: 

Sxi = ucos(Q • x) + vcos(Q • x) (19) 

where Q is the propagation vector (pitch) of the spiral, 
and corresponds to a minimum of J™'". In spin space, 
the plane of the spiral is fixed by two orthonormal vec- 
tors u and V. When m = 1, it is only when J™'" admits 
several degenerate minima in the Brillouin zone that ad- 
ditional non-spiral (and possibly non-planar) GS may be 
constructed. If the lattice has more than one site per unit 
cell, an attempt to construct a spiral with a pitch corre- 
sponding to the smallest eigenvalue J(Q)'"'" will gener- 
ally not lead to a physical spin configuration with fixed 
spin length Sf^ = 1 at every site. We will however see 
in the next section that for some models, a non-planar 
regular state may reach the lower bound, whereas all the 
spiral states are energetically higher. 

C. Variational phase diagrams of Heisenberg 
models on the kagome and hexagonal lattices 

In this section we comment the phase diagrams of Ji- 
J2-Jz{-J''i) Heisenberg models on the kagome and hexag- 
onal lattices. J„ is the interaction between n*'* neigh- 
bors. On the kagome lattice, there are two types of third 
neighbors depicted in Fig. 8(a), and thus two coupling 
constants J3 and J3. 

For each set of parameters, we determined the regular 
state with the lowest energy, the SS of Sec. IV D with 
the lowest energy, and the value of the lower bound on 
the energy (the energies of regular states are given in 
Figs. 4 and 5). The results on these two lattices are 
described in Figs. 7 and 8. Such phase diagrams are a 
priori variational. However, it turns out than in all the 
colored (white included, grey and black excluded) regions 
of Figs. 7 and 8, the regular state with the lowest energy 
reaches the rigorous lower energy bound of Eq. (18). This 
demonstrates that (at least) one GS is regular in these 
regions of the parameter space. In the grey areas, the 
energy lower bound is not reached, but the regular near- 
by state could be a GS as no SS has a lower energy. 
In the black areas, the GS is not regular: some SS is 
energetically lower (but sometimes still higher than the 
lower bound). All regular states (excepted those from 
continua) appear in some area of the presented phase dia- 
grams. This shows that these states are good candidates 
as variational GSs. The absence of regular states of a 
continuum in an Heisenberg model is easily understood. 
The energy E of any regular state c belonging to a con- 
tinuum cannot be lower than the energies Ei and E2 of 
the two states ci and C2 between which it interpolates. 
One (at least) of the two states, say ci, is colinear along 
a direction n. The C2 spins are then perpendicular to 
n. Let 9 be the angle between the spins of the contin- 
uum state and n. Then Si = S^^ cos 6 + S^^ sin and 
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FIG. 7: (Color online) Phase diagram of the J1-J2-J3 Heisen- 
berg model on the honeycomb lattice. Labels refer to the reg- 
ular states described in Fig. 5. In each colored region (black 
excluded), the regular state is an exact GS. In the black re- 
gion a generalized spiral state (SS) has an energy stricly lower 
than the regular states, but the actual GS energy might still 
be lower. 




(a) Definition of the coupling constants of the model. 
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FIG. 8: (Color online) Phase diagram of the J1-J2-J3-J3 
model on the kagome lattice. In each colored region (white 
included, grey and black excluded), the regular state is an ex- 
act GS. Labels refer to the regular states described in Fig. 4. 
In the grey regions, the near-by regular state does not reach 
the lower bound of Sec. VI B but no SS state is energetically 
lower. In the black regions, a SS has a lower energy than the 
regular states, but the actual GS might yet be lower. 



the energy reads E = E2 + {Ei — E2) cos^ 9. Thus, E 
is in between Ei and E2 and is never strictly the lowest 
energy." 

We will now address the possible degeneracies of reg- 
ular tridimensionnal spin states in these models. On the 
hexagonal lattice, our phase diagram is in agreement with 
Ref. One should nevertheless notice that the regular 
tridimensionnal orders (tetrahedral and cubic states) are 
degenerate with colinear non regular states. These last 
states have a higher density of soft excitations (larger 
energy wells in the phase space landscape) and will al- 
ways win as soon as (thermal or quantum) fluctuations 
are introduced (order by disorder mechanism^ '"' '). How- 
ever the non planar configurations could be stabilized by 
quartic or ring-exchange interactions. 

On the kagome lattice (Fig. 8) the occurrence of the 
cuboc2 (Fig. 4(f)) for J1-J2 interactions and of the 
cubocl (Fig. 4(e)) for J1-J3 interactions has already 
been reported. These two states are not degenerated with 
SS and are to our knowledge unique and stable GS of the 
model. To our knowledge, the octahedral state has not 
been found before, but this state has the same energy as 
a continuum of non SS states including colinear states, 
and it will be destabilized by any fluctuation. 



D. Square and triangular lattices: Phase diagrams 
of Heisenberg versus ring-exchange models 

In this section we will comment the phase diagram 
of the Heisenberg models (Eq. (1)) on the square and 
triangular lattices and display the effect of 4-spin ring- 
exchange {J1-J2-K) on these two lattices. The J1-J2-K 
model is deflned as: 



i? = ^J(|x,-x,|)S,.S, 
+(S, • Si){Sj • Sfe) - ( 

+ Sj • Sfc -I- Sfc • S; -I- S 



K 



^ ((s, •s,)(Sfc.so 

>z ■ Sfc)(Sj • Si) -f Si • Sj 
•S, + S, -Sfe + S, -S,) (20) 



where the sum in the K term runs on rhombi i^ji k,L' 
This model encompasses flrst and second neighbor Ji and 
J2 couplings and a K ring-exchange term which intro- 
duces quartic interactions as well as modiflcations of first 
and second neighbor Heisenberg interactions. The phase 
diagramms are displayed in Figs: 9 and 10. 

In the J1-J2-J3 Heisenberg phase diagrams on the 
square and triangular lattice, all regular states that do 
not belong to continua do appear as exact GSs in some 
parts of the phase diagrams (colored regions - black ex- 
cepted of Figs. 9(a), 9(b), 10(a) and 10(b)). In black 
regions, SS are more stable than regular states. As these 
lattices are Bravais lattices, we know how to reach the 
lower bond of Sec. VI B thanks to a spiral state. The 
orthogonal state on the square lattice and the tetrahe- 
dral state on the triangular lattice (Fig. 3(b) and 6(c)) 
are degenerate with spiral states including colinear states 
with 2 spins (up, down) in the magnetic unit cell, which 
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FIG. 9: (Color online) Phase diagrams on the square lattice 
with J1-J2-J3 Heisenberg interactions (top line) and J1-J2-K 
model (bottom line). Labels refer to regular states defined in 
Fig. 6. In each colored region (black excepted), the regular 
state has the lowest energy of the set of all regular and SS 
states. In the black regions, a SS has a lower energy than 
the regular states. For pure Heisenberg interactions, we know 
that we obtain the GS energy, but for non Heisenberg interac- 
tions the actual GS might be lower. In the J1-J2-J3 model the 
coplanar (orthogonal 4-sublattice) phase is degenerate with 
non regular colinear states, which will win upon introduc- 
tions of fluctuations. A contrario the coplanar (orthogonal 
4-sublattice) phase is stable in a large range of parameters in 
the J1-J2-K model. 



will win upon introduction of fluctuations. On a large 
part of the phase diagram on the square lattice (spirals 
excepted) the spins are thus colinear. 

The presence of a 4-spin ring exchange on the square 
lattice gives richer phase diagramms (Figs. 9(c) and 9(d)) 
with the appearance of states from continua. We re- 
call that these phase diagrams are variational and give 
the minimal energy state among the regular and the SS 
states. A dominant 4-spin ring exchange stabilizes the 
orthogonal 4-sublattice coplanar antiferromagnet , which 
is known to be robust to large quantum fluctuations. 
One of these phases belongs to a continuum: the V states 
(Fig: 6(d)). Part of this phase diagram on the square lat- 
tice has been known for a long time for the Ji-K model,-" 
but the effect of a second neighbor interaction leads to 
new phases that might be interesting in various respects. 

The J1-J2-K phase diagramm on the triangular lattice 
(Figs. 10(c) and 10(d)) exhibits all the regular phases 
that can be constructed on this lattice. In that model, 
large ring-exchange stabilizes the tetrahedral chiral phase 



FIG. 10: (Golor online) Phase diagrams on the triangular 
lattice with J1-J2-J3 Heisenberg interactions (top line) and 
J1-J2-K model (bottom line). Labels refer to regular states 
defined in Fig. 3. In each colored region (black excepted), 
the regular state has the lowest energy of the set of all regular 
and SS states. In the black regions, a SS has a lower energy 
than the regular states. For pure Heisenberg interactions, we 
know that we obtain the GS energy, but for non Heisenberg 
interactions the actual GS might be lower. 



studied by Momoi and co-workers."''^ The presence of 
large parts of the phase diagramms with planar or 3- 
dimensional order parameter at T = 0, and of points 
where a large number of classical phases are in compe- 
tition, could give interesting hints in the quest of exotic 
quantum phases. 



E. Finite temperature phase transitions in 
two-dimensions 



In two-dimensions, the Mermin-Wagner" ' theorem in- 
sures that continuous symmetries cannot be sponta- 
neously broken at flnite temperature. It does however 
not prevent discrete symmetries to be broken. Indeed, 
some flnite temperature phase transitions associated to 
discrete symmetries have been found in classical 0(3) 
models: lattice symmetry breaking in the J1-J2 and J1-J3 
models on the square lattice, " chiral symmetry break- 
ing in a ring-exchange model on the triangular lattice^ ''"^ 
and in a J1-J2 model on the kagome lattice. ^'''^'^ 

What should be expected in a system where the GS is 
a regular state c ? Let us flrst consider the case where c is 
not chiral, that is when the spin inversion S — >■ — S gives 
a state c' which can also be obtained from c by a rotation 
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in 5*0(3). At an infinitesimally small temperature, the 
rotational symmetry is restored and the statistical ensem- 
ble is that of all the (regular) states obtained from c by 
SO{3) rotations. The thermal average of an observable 
is therefore also an average over iS'0(3) rotations. Now, 
if we compare an observable O and the same observable 
after a lattice symmetry X, we will get the same aver- 
age (for regular states, the effect of X can be absorbed 
by a rotation). So not only the rotational symmetries, 
but all the lattice symmetries are restored at T = 0"''. 
The simplest scenario is therefore a complete absence of 
symmetry-breaking phase transition from T = 0^ up to 
T = oo. Now, for a chiral state, the thermal fluctua- 
tions will only partially restore the 0(3) symmetry of 
the model, and a chiral phase transition - possibly ac- 
companied by some lattice symmetry breaking - should 
be expected. From this point of view, a classical sys- 
tem in two-dimensions with no finite-temperature phase 
transition is likely to have a regular and non-chiral GS. 



VII. CONCLUSION 

Based on symmetry considerations (and on an analogy 
with Wen's ' classification of quantum spin liquids using 
the concept of PSG), we introduced a family of classical 
(antiferro-) magnetic structures, dubbed "regular" states. 
They can be constructed in a systematic way for any 
lattice, based on the method explained in Sec. III. We 
found that these states are often good variational states 
to study the zero-temperature phase diagram of "com- 
plex" problems (non-Bravais lattice and/or multiple spin 
interactions for instance). In many cases, one of the regu- 
lar state is found to reach a lower energy bound, allowing 
to show that it is a GS. 

We note that, although one can always find a planar 
GS in Heisenberg models on a Bravais lattice, non-planar 
antiferromagnetic spin structures with many sublattices 
are rather common in presence of competing interactions, 
non quadratic spin interactions and non-Bravais lattices. 
As mentioned in the introduction, we believe this ap- 
proach may find an application in the study of real mag- 
netic compounds where the (equal time) spin-spin corre- 
lations are measured, but the strength and range of the 
magnetic exchange interactions are not known. 

We have studied the case where the spin manifold 
^ = is that of a three-component spin (unit vector), 
but other manifolds could be investigated using the same 
approach. For instance, nematic regular orders would be 
obtained with A = and Ss — S0{3). 
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Appendix A: Powder-averaged structure factors of 
regular states 

Equal time spin-spin correlations partially characterize 
a spin state and are independent of the energetic prop- 
erties of the system. Equal time structure factors can 
thus be analytically calculated on regular states to form 
a set of reference neutron scattering results. They can be 
used to analyze measurements done on compounds with 
unknown GS. We define the equal time structure factor 
5(Q) of a state as 

S'(Q) oc J2 e~''^^^'"''^^S, • Sj, (Al) 

where x.^ is the position vector of the site i. The pro- 
portionality factor is adjusted to verify the sum rule 
5'(Q) = 1. For perfect long-range orders, S{Q) is 
zero everywhere except for a finite number of Q where 
Bragg peaks are present. They are broadened when 
chemical defects, non zero temperature or quantum fluc- 
tuations are taken into account. 

When only powders are realisable, one can measure the 
powder equal time structure factor S'dQj). It is the aver- 
age of 5(|Q| sin6'(u cos -I- V sin ■)/))) over all the possible 
3d orientations of Q, where 9, tp are the spherical coor- 
dinates angles of Q in the orthonormal basis (u, v, uA v) 
with u, V in the sample plane. Thus 

^(IQI) oc / d\ Q(IQ|-Iql) 5(q), (A2) 

where 9 is the Heaviside step function and q browses the 
reciprocal 2d space. 

The equal time structure factors 5(Q) were given in 
Fig. 3, 4, 5 and 6 for the regular states on the triangular, 
kagome, honeycomb and square lattices. The powder- 
averaged equal time structure factors S'dQj) on the tri- 
angular and kagome lattices are shown in Fig. 11 and 12. 



Appendix B: Analogy with Wen's Projective 
symmetry groups (quantum spin models) 

For quantum spin-^ Heisenberg models, a standard 
mean-fleld approximation consists in expressing the spin 
operators in term of fermionic operators fi^, where i is 
a lattice site and a is the spin ±1/2. A mean-fleld 

decoupling based on some bond parameters r]ij and 
(notations and details to be found in in Ref. M) can then 
be performed to make the Hamiltonian quadratic in the 
fermionic operators. 
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FIG. 11: Powder-averaged equal time structure factors 5(1 Ql) 
of the regular states on the triangular lattice (|Q| is in units 
of 27r, S(IQj) in arbitrary units). 
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(d) Octahedral state (e) Cubocl state (f) Cuboc2 state 

FIG. 12: Powder-averaged equal time structure factors 5(|Q|) 
of the regular states on the kagome lattice (|Q| is in units of 
27r, 5(|Q|) in arbitrary units). 



This theory has a local SU{2) gauge invariance. The 



set of gauge transformations is denoted by $. Physi- 
cal quantities, which can be expressed using spin opera- 
tors, are unaffected by a gauge transformation, although 
and are generally modified. A mean-field state is 
characterized by a set of rjij and S^ij values, called Ansatz. 
Two mean-field states do have the same physical observ- 
ables if they are related by a gauge transformation. The 
group of transformations (lattice, gauge and combined 
transformations) that do not modify an Ansatz is called 
the projective symmetry group (PSG). Its subgroup of 
pure gauge transformations is called the invariance gauge 
group (IGG). 

One may be interested in states for which all the physi- 
cal quantities are invariant under the lattice symmetries. 
To classify these "uniform" states, one can first fix the 
IGG and then look for the "algebraic" PSG which obey 
the constraints derived from the algebraic structure of 
lattice symmetry group Sl ■ The actual Ansatze can then 
be constructed. 

Clearly, there is a close correspondence between the 
construction of regular states discussed in this paper, and 
that of symmetric Ansatze. This correspondence is sum- 
marized in Tab. I. 





Classical spin 
models 


Quantum 
Mean-field 


State 


Regular state 


Physically 
symmetric Ansatz 


Internal symmetry 
group 


Ss (global spin 
rotation, etc.) 


$ (local gauge 
transformations) 


Symmetry group 
of a state 




PSG 


Unbroken internal 
symmetries 




IGG 



TABLE I: Analogy between the construction of regular states 
and that of symmetric Ansatze in :. 
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Once ah the regular states have been constructed for given 

lattice and spin symmetries (using a simple group theoretic 

construction, as explained in Sec. Ill), one can directly 

compare the energies they have for a given microscopic 

Hamiltonian. 

A case where Sh Ss x Sl is the antiferromagnetic 
square lattice with a site-dependent magnetic field taking 
two opposite values on each sublattice. The spin inversion 
Si — >■ Si is not in Ss, the translation by one lattice spac- 
ing is not in Sl, but the composition of both is in Sh- The 



theory developed in this paper can however be used in this 
case by replacing Sl by Sh/Ss- 

Notice that nearest neighbors on the lattice do not neces- 
sarly correspond to nearest neighbor spin directions in spin 
space. 

Again, the plaquettes of the lattice need not to map to the 
faces of the polyhedron. 

This relation is particularly easy to visualize in the case of 
the tetrahedral state on the triangular lattice, since both 
the lattice and the polydedron E have triangular plaque- 
ttes/faces: one can put a tetrahedron with a face posed 
onto a lattice face. Then, one roll the tetrahedron over the 
lattice to obtain a spin direction at each lattice site. Notice 
that such a construction would not work with a cube on the 
square lattice (and indeed, there is no such eight-sublattice 
regular state on the square lattice, see Sec. IV C). 
In the presence of an external magnetic field h and if 
a one-dimensional representation included in G is ferro- 
magnetic (as it is the case for some umbrella's and for 
the V-states), n aligns on h. The energy then reads E = 
E2 + {El — E2) cos^ 6 — hcosS and an umbrella state be- 
comes stationnary. It is well known that such structure can 
be the GSs in presence of a magnetic field. 



